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ABSTRACT 

The 14 C production of shock-accelerated particles is calculated in terms of 
the total energy released in energetic particles. The recently reported 1.2% jump 
in the 14 C content of the atmosphere in the year C.E. 775, it is found, would 
require > 10 34 erg in energetic particles, less than first estimates but far more 
than any known solar flare on record. It is noted that the superflare from a 
large comet (comparable to C/Hale-Bopp) colliding with the sun could produce 
shock-accelerated GeV cosmic rays in the solar corona and/or solar wind, and 
possibly account for the CE 775 event. Several additional predictions of cometary 
encounters with the sun and other stars may be observable in the future. 

Subject headings: Comets; flares 
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1. Introduction 

Comets, asteroids and giant solar flares each pose dangers. If a comet, or its coma or 
tail were to sufficiently rattle the Earth's magnetosphere, the electromagnetic disturbance so 
induced could threaten modern civilization, which depends on functioning microelectronics. 
A sufficiently large sun-grazing comet (or asteroid) R > 10 7 cm, would have a mass of 
10 2L5 g and contain over 10 36 ergs, more than two orders of magnitude more kinetic energy 
near the sun than a reasonable estimate for the energy in the Carrington solar flare of 1859 
(which damaged early telegraph lines), and its energy release near the sun's surface could 
traumatize the Earth with UV exposure and electromagnetic disturbance. A mid-size to 
large comet (R > 3 km) impacting Earth could deposit > 10 30 ergs into the ocean, enough 
to supersaturate the Earth's atmosphere with water vapor, leading to something resembling 
legendary floods, and could dramatically heat at least parts of the atmosphere. 

Sun-grazing comets are a fact of life (e.g. Schrijver, et. al 2012, Sekanina and Chodas, 
2012). The possibility clearly exists that some of them could be, at some stage, quite 
large (M 3> 10 19 g) and collide with the sun, causing an explosive release of more than 
10 34 erg in energy (Brown et al., 2011). The lack of any known, reliable record of a major 
cataclysm associated with such past events could be interpreted - depending on the size and 
kinetic energy of the comet - as evidence that they were less conspicuous to a pre-electronic 
civilization than one might suppose, or that superstitions interpreting comets as bad omens 
had some historical basis. In any case, they could be far more disruptive to post-20th 
century civilization, and their event rate, even if small, is worthy of study. 

Here we consider whether a) a giant solar flare, or b) the close approach of a large 
comet to the sun could have occurred in the year 775, when the levels of 14 C rose by 1.2 
percent within a year or so (Miyake, et al. 2012). This rise would require 10 years worth of 
normal cosmic ray exposure within one year. While it could have been due to an extremely 
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large but otherwise normal solar flare, its statistical deviation from other 14 C rises on record 
motivates us to consider whether a different type of event from normal solar flares could 
produce the 14 C enhancement. If it was indeed due to a giant comet closely encountering 
the sun, it would be evidence that such an event is survivable, and quantitative estimates 
of the energy it released are desirable. By the same token, it is worth considering whether 
a less close approach could have produced the enhancement by tapping the energy in the 
solar wind to produce a temporary rise in energetic particles. 



2. Particle Acceleration and 14 C Production 

The energetic particle spectrum F(p)dlnp at momentum p for shock-accelerated ions 
at the shock is given as 
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(Eichler, 1981a). Here r is the compression ratio of the shock, 1//2 = pvc/2 1 ^ 2 nr]ZeB 

is the geometric mean of the parallel and perpendicular diffusion coefficients of ions of 
momentum p, velocity v, and charge Ze in the magnetic field B, R s is the radius of the shock, 
and u s is its velocity relative to the pre-shock fluid. The spectrum cuts off exponentially at 
energy E > E Q = i]ZeBR s u s /c, where the dimensionless coefficient r\ is about 1/3, based 
on observations of the energetic particles at the Earth's bow shock (Ellison and Mobius, 
1987, Chang et al, 2001) at which R s ~ 6 • 10 9 cm and E Q ~ 36Z KeV. For pv < E a , the 
integral spectral index s is -1 for a strong shock of compression ratio 4 in the test particle 
approximation. 

The spectrum of escaping particles, in the approximation of steady state, is proportional 
to their rate of production at the shock, which is proportional to u s (l — l/r)(dF/dlnp) /3. 
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This expression includes particle escape by both convection and diffusion to a free streaming 
boundary. We assume that the particles mostly responsible for the 14 C, i.e. the most 
energetic, freely stream from a sunward acceleration site towards Earth and precipitate 
onto its polar caps. Particles trapped in the expanding flow are those at low energy, and 
their adiabatic losses can be recycled into the acceleration of the expanding blast wave. 

The 14 C production rate per unit particle energy J Ejj^dlnp in the Earth's atmosphere 
is independent of the normalization of F and is given by 



where E is the kinetic energy, hereafter expressed in units of m p c 2 , Y(E)/ir, the neutron 
yield per primary cosmic ray of energy E, can be approximated as 4E 2 - 35 / (E 2 + E - 35 ) 
(Kovaltsov et a!., 2012 ), Q,(E)/(2tt) ~ 0.13E 1 / 4 is the solid angle of the magnetic polar 
cap whose geomagnetic cutoff is E. (Here we have assumed that the angular distribution of 
the energetic particles is energy independent.) Using equations (T25 and (j2J), we numerically 
calculate the total 14 C yield per unit energy as a function of the parameter E Q . For 
reference, note that the most energy-efficient ion energy for making 14 C is at the maximum 
of Q(E)Y(E)/E, and that the yield is ~ 0.3 14 C atoms per m p c 2 of kinetic energy. The 
results for an actual shock-accelerated spectrum over a wide range of E a are plotted in 
figure 1, showing that over a wide range of E a , the 14 C production efficiency e is within a 
factor of several of the maximum. 

An episode in which 6 • 10 8 14 C atoms/cm 2 are produced at Earth then requires a 
(kinetic) energy fluence at Earth of ~ 6e _1 • 10 8 m p c 2 /cm 2 . A giant particle acceleration 
event near the sun that spewed out energetic particles over a solid angle of ~ 2n would thus 
have needed to produce 8.4 • 10 35 e _1 m p c 2 ~ 1.3 • 10 33 e _1 ergs in energetic particles. The total 
energy of the event, of course, could be more, though several observations in the heliosphere 
(e.g. Eichler, 1981b, Ellison & Mobius, 1987, Decker et. al, 2008) support the theoretical 
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assertion (Eichler, 1979, Eichler, 1985, Ellison & Eichler, 1985) that energetic particles can 
contain a significant fraction of the total energy in a collisionless blast. The value of the 
cutoff energy E Q = eBRu s /3c obtained by using solar parameters - B > 0.3 G, R ~ 6 • 10 10 
cm, u s /c > 1CT 3 - is comfortably above 1 GeV, so, even if the lateral extent of the shock 
is somewhat less than the above value of R, we may reasonably assume a value of e in the 
range 0.1 to 0.2. A superflare yielding ~ 10 34 ergs in energetic particles could thus account 
for the 14 C enhancement of 775. For E Q > 0.3m p c 2 , this energy requirement is less than the 
estimate in Miyake et al. (2012), presumably because the particle spectrum predicted here 
is more efficient for 14 C production than the hardest recorded solar flare spectra, which 
was adopted by those authors. Nevertheless, our minimum energy estimate is far enough 
beyond the energy of recorded flares that we are motivated to consider whether the event 
could have been of cometary origin. 

A superflare of energy 10 34 erg could be produced by a solar encounter of a comet (or 
asteroid) of mass > 10 19 g, which is nearly the mass of C/Hale-Bopp (Weissman, 2007). It 
has been estimated that comet Hale-Bopp has a 0.15 chance of eventually colliding with 
the sun within several hundred orbits (Bailey, Chanbers, & Han, 1992, Bailey et al, 1996). 
Given that its most recent apparition was less than 2 decades ago, let us assume that a 
comet this large appears about twice per century. Multiplying by 0.15, and assuming each 
one reappears ~ 300 times before colliding with the sun we roughly estimate, on the basis 
of this very small number statistic, a minimum solar collision rate of > 10~ 5 yr -1 . 

The above estimate is a lower limit based on a single, long period comet. Of the 
~ 45, 000 Centaurs estimated to exist, whose orbital stability times are of order 10 6 to 10 7 
yr., many are thought to scatter off Jupiter or outer solar system bodies into the inner solar 
system, where some meet their demise by crashing into the sun (Horner, Evans & Bailey, 
2004). So a new Centaur enters the inner solar system every ~ 10 2 years or so. If ~ 0.06 
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of them collide with the sun (Levison & Duncan, 1994), it would be plausible to expect 
such an event every ~ 1.5 millenia, comparable to the time scale spanned by the various 
tree-ring and ice core data sets. 

The Kepler telescope data set reported superflares from 365 stars that displayed 
brightness variations (> 10~ 3 ) indicative of a magnetic origin (Maehara et al, 2012). In 
addition, there were 9 superflares from stars, out of a sample of 4.5 • 10 4 stars that showed 
otherwise steady emission to within instrumental error (Maehara, private communication). 
The upper limit on the event rate for stars without periodic luminosity fluctuations greater 
than 10~ 3 is thus about 3 • 10~ 4 yr^ 1 for all stars, and the question is still open as to 
what fraction have planetary systems capable of driving comets into their host stars. We 
conclude that the hypothesis of a C/Hale-Bopp -size comet hitting the sun every several 
millenia is consistent with present observations. 



3. Cometary Bow Shocks 

We have also considered a scenario in which the comet converts solar wind energy to 
energetic particle energy at the bow shock made by its coma. Because the radial B field 
decreases with distance D from the sun as B(D) ~ 10~ 5 5 5 [1A.U./D] 2 G, b 5 ~ 1, while the 
radius R s of the bow shock could increase at most as D, it is clear that E Q decreases with 
D, and can never be more than r/eBDu sw /c ~ 60[lA.U./D]MeV. For E Q to be above 300 
MeV, D must be at most 0.2 A.U., and the amount of time that the comet would spend 
within this distance can easily be shown to be about 5t — 4 • 10 5 s. Assuming the freely 
streaming energetic particles and the solar wind suffer the same inverse square dilution in 
getting from the sunward comet to Earth, it suffices to consider the kinetic energy flux of 
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the solar win«j^ (sw) at Earth, p sw u 3 s /2, times St. This sets a maximum energy fluence in 
energetic particles at Earth of p sw u^5t/2 = 7.2 • 10 7 (n slu /5cm 3 )(u s /4 ■ 10 7 cms _1 ) 3 m p c 2 /cm 2 
and a maximum 14 C production of < 10 7 cm -2 . This upper limit falls short of the inferred 
14 C production by more than an order of magnitude. A train of N large comet fragments 
(N > 10 2 ), all from a single progenitor that fragmented, could each play out the scenario 
N times within a year, and thus enhance the above estimate. Moreover, the pitch angle 
of escaped particles arriving from a sunward source could be biased toward the parallel 
direction, thus decreasing the amount of mirroring at the poles relative to Galactic cosmic 
rays, and increasing the fraction of precipitating particles. Nevertheless, we believe the 
original scenario - a C/ Hale- Bopp- size comet crashing into the sun - to be the more 
conservative, plausible scenario. 



4. Further Discussion 

If the most likely explanation of the CE. 775 event is a superflare at the sun's surface 
- as opposed to an event within the solar wind per se - then the lack of any record of 
devastation of any sort at that time is reassuring, though not in regard to safety of power 
grids and satellites. On the other hand, it implies a huge blast and energetic particle 
flux, delivered impulsively, at Earth. If it were to repeat in the modern era, it could have 
devastating consequences for power grids and satellite electronics. 

Most sun-grazing comets are of low mass, M < 10 12 g, and carry < 2 • 10 27 erg of 
kinetic energy into the sun. This is probably not enough to be detectable in high energy 
particles or their secondaries. Comet Lovejoy, the one known sun-grazer that was large 
enough to survive the perihelion of a close solar encounter, probably had a surface area of 



with a modest correction for the comet's motion, which adds to u s 
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about 10 10 cm 2 and deposited 10 13 g to 10 14 g (Sekanina and Chodas, 2012) and 2 • 10 28 to 
2 • 10 29 ergs during its passage through the corona. Assuming a coma size of 10 9 i? 9 cm, a 
magnetic field of B G, and a shock velocity of 10 _2 ' 5 /3_ 2 .5C, the value of the cutoff energy 
E Q is E Q = 300i? 9 -Bo/3-2.5 MeV. It is therefore conceivable that high energy ions, and/or 
secondary neutrons and gamma rays were produced at detectable levels during the passage. 
Results from the IMPACT mission of STEREO are therefore awaited at the time of this 
writing. It is possible that at some time in the near future a sun-grazing comet will produce 
an energetic particle event that will teach us a great deal. We have also noted that a more 
extensive data set from the Kepler Observatory could reveal non-magnetic superflares on 
solar-type stars. 

Finally, the 10 Be enhancement during the same event sets an independent constraint 
on the spectrum of energetic particles that caused the CE 775 event (Miyake et al. 2012). 
Though beyond the scope of this paper, a careful analysis of the 10 Be data during that time 
could test our rough estimate for the spectral cut-off parameter E Q > 1/3 GeV, and the 
attendant value of the 14 C production efficiency e. 
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Dr. M. Eichler, for helpful conversations, and the support of the Israel-U.S. Binational 
Science Foundation, the Israeli Science Foundation, and the Joan and Robert Arnow Chair 
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